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ABSTRACT 

Transverse oscillations of coronal loops are often observed and have been theo- 
retically interpreted as kink magnetohydrodynamic (MHD) modes. Numerical simu- 
lations by Terradas et al. (2008, ApJ 687, LI 15) suggest that shear flows generated 
at the loop boundary during kink oscillations could give rise to a Kelvin-Helmholtz 
instability (KHI). Here, we investigate the linear stage of the KHI in a cylindrical 
magnetic flux tube in the presence of azimuthal shear motions. We consider the basic, 
linearized MHD equations in the /3 = approximation, and apply them to a straight 
and homogeneous cylindrical flux tube model embedded in a coronal environment. 
Azimuthal shear flows with a sharp jump of the velocity at the cylinder boundary are 
included in the model. We obtain an analytical expression for the dispersion relation 
of the unstable MHD modes supported by the configuration, and compute analytical 
approximations of the critical velocity shear and the KHI growth rate in the thin tube 
limit. A parametric study of the KHI growth rates is performed by numerically solv- 
ing the full dispersion relation. We find that fluting-like modes can develop a KHI in 
time-scales comparable to the period of kink oscillations of the flux tube. The KHI 
growth rates increase with the value of the azimuthal wavenumber and decrease with 
the longitudinal wavenumber. However, the presence of a small azimuthal component 
of the magnetic field can suppress the KHI. Azimuthal motions related to kink oscil- 
lations of untwisted coronal loops may trigger a KHI, but this phenomenon has not 
been observed to date. We propose that the azimuthal component of the magnetic field 
is responsible for suppressing the KHI in a stable coronal loop. The required twist is 
small enough to prevent the development of the pinch instability. 
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1. INTRODUCTION 



The phenomenon of transverse coronal loop oscil lations has received much attention during 
the la st decade since the first observational reports (e.g. jAschwanden et al. iNakariakov et alj 



and th eir subsequent theoretical interpretation in terms of kink rnagnetohydrodynamic ( MHD) 
modes (e.g.. INakariakov & OfmanlboOll : iRuderman & Robertsll2002l : looossens et aDl2002l) . The 
reader is referred to recent reviews, and refere nces therein, regarding the m odeling of transverse 
loop oscill ations in terms of MHD eigenmodes (IRuderman & Erdelvill2009h . their time-dependent 
excitation ( Terradas 2009 ). and their non-linear evolution ( Ofmanl2009 ). 



Terradas et al.l (|2008bh studied 3D numerical simulations of non-linear kink oscillations of a 
straight, untwisted coronal magnetic flux tube. These authors found that, after kink motions are 
excited, small length-scale disturbances appear at the tube boundary and grow rapidly in time. The 
inclusion of a transversely inhomogeneous transitional layer between the flux tube and the exter- 
nal medium does not suppress the apparition of these small-scal e disturbances, bu t they develop 
more slowly than in the sharp transition situation. Although Terradas et al. ( 20081 ^) qualitative in- 
terpreted this phenomenon in terms of a Kelvin-Helmholtz instability (KHI), which would be trig- 
gered by shear flows at the tube boundary, a more detailed investigation of the instability regimes 
and growth rates is needed. This can be performed in two diff"erent ways. The first option is to study 
the n on-linear evolution of t he KHI by numerically solving the full 3D non-linear MHD equations 



as m iTerradas et al.l (l2008b() . but with a very large spatial resolution able to correctly describe the 
small spatial-scales generated at the tube boundary. This procedure requires an enormous compu- 
tational effort. The second option is to attack the problem analytically by restricting ourselves to 
the initial stage of the KHI and considering its linear evolution in a simplified model that includes 
the basic ingredients of the actual situation. By means of an analytical work, a more in-depth 
investigation of the effects and parameters involved in the generation of the KHI is possible. We 
adopt this second approach in the present work. 

To do so, we have to identify the basic ingredients that the model must contain. We consider 
a straight cylindrical flux tube embedded in a coronal environment to represent a coronal loop. It 
is well-known that for kink oscillations of the flux tube, the azimuthal componen t of the velocity 



pertu rbation has a finite jump at the cylinder boundary (see, e.g.. Figure 1(b) of iGoossens et al. 



20091). We claim t hat this jump of the az imuthal velocity is responsible for the observed KHI in 
the simulations by iTerradas et al.l (l2008bl) . So, we take the presence of azimuthal shear flows at 
the cylinder into account to represent those azimuthal motions generated by the kink oscillations 
themselves. However, to consider the actual azimuthal velocity profile of the kink mode perturba- 
tion (expressed in terms of Bessel functions) makes the problem not analytically tractable. Instead, 
we assume here a simpler profile for the azimuthal velocity which allows us to proceed analyti- 
cally, but keeping the presence of a sharp velocity discontinuity at the cylinder boundary as the 
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key ingredient of the real situation. We also neglect the dependence of the flow on the longitudinal 
direction. In addition, we assume that the azimuthal shear flows are not time-dependent. This 
condition restricts ourselves to KHI growth rates much smaller than the kink mode period, but it 
allows us to perform an analysis in terms of normal modes. 

The effect of mass flow on coronal loop oscillations h as been studied by a number of au 



thors , but they mainly focused on longitudinal flow s (e.g., iGoossens et al.l 1 19921 : Erdelyi et al. 



19951 : iTerra-Homem et al.ll2003l : iTerradas et al.ll2009l. among other works) . Applications to fila 
ment thread oscillations have also been performed dSoler et al.ll2008L |2009|) . Concerning the sta 
bility of these flows, both the KHI and the resonant instability have been inves t igated (e.g.. 



1983 



Andries et al.l l2000l : lAndries & GoossensI 1200 ll : Erdelyi & TaroyanI l2003l : iHolzwarth et al 



Rae 



20071). Regarding azimuthal shear flows, there are several previous investigations in the context of 
astrophysical and laboratory plasmas that are relevant for the present application. A very compre- 
hensive investig ation of a similar config ur ation to that adopted here but in the context of tokamaks 
can be found in iBondeson et al.l (Il987|) . iBodo et al.l (|l989l 119961) studied the stability of cylin- 
drical jets in the presence of both azi muthal and longitu dinal flows. While they first restricted 



themselves to the axisymmetric case in bodo et all Jl989 t). they later extended their investigation 



to non- axisymmetric perturbations in IBodo et al.l (Il996|) . These authors obtained an analytical ex- 
pression for the dispersion relation of the unstable modes an d computed the instability growth rates 
under astrophysical jet conditions. l Ogilvie&Pring leld 19961) studied the s tability of a differentially 
rotating accretion disk in the presence of an azimuthal magnetic field, and lKolesnikov et al.l (|2004r) 
investigated the KHI due to helical flow around a photospheric magnetic flux tube. To our knowl- 
edge, no similar analytical work applied to the context of kink oscillations of coronal loops is found 
in the literature. 



2. BASIC EQUATIONS 

Our configuration is composed of a homogeneous, straight, and cylindrical magnetic flux tube 
of radius a embedded in a homogeneous coronal environment. The geometry of the model allows 
us to use cylindrical coordinates, namely r, ip, and z, for the radial, azimuthal, and longitudinal co- 
ordinates, respectively. In the following expressions, a subscript indicates equilibrium quantities 
while we use subscripts i and e to explicitly denote internal and external quantities, respectively. 
The magnetic field is taken homogeneous and orientated along the cylinder axis, Bq = Bqc^, with 
Bq constant everywhere. We assume a step profile for the density, with pi and pe the internal and 
external densities, respectively. We consider an azimuthal shear flow in the equilibrium configura- 
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tion, Uo = Uqc^, with Uq given by 



0, if r > a. 



Thus, voa is the velocity jump at the cylinder boundary. We are aware that the selected flow profile 
is substantially different from the kink mode azimuthal velocity. However, our flow profile keeps 
the key ingredient of the actual kink mode profile, i.e., the presence of a sharp jump at the cylinder 
edge. The KHI is mainly governed by this jump of the azimuthal velocity, while the precise flow 
profile within the flux tube would only have a minor influence on the growth rate s . In a ddition, 
the present profile allows us to proceed analytically. Here, and as in iBodo et all d 19891 . E996), 



we assume that v^a «; Vai, with vai the internal Alfven speed (defined below), which allows 
us to neglect the eff"ect of the centrifugal force over the equilibrium total pressure gradient and 
therefore the equilibrium quantities can be taken constant in the radi al direction. Since part of the 



following treatment is equivalent to that of iBodo et al.l (|1989l Il996|) . we will refer to their works 
when appropriate. 

We consider the linearized MHD equations for an ideal plasma in the /? = approximation, 
namely 

Po + Uo ■ Vvi + vi ■ VUoj = (V X Bi) X Bq, (2) 

^ = Vx(UoxBi) + Vx(vixBo), (3) 

along with the condition V • Bi = 0. In these equations Bi = {Br, B^, B^, Vi = (v^, v^, v-j are the 
components of the magnetic field and velocity perturbations, respectively, while fi is the vacuum 
magnetic permeability. Next, we write the perturbations proportional to exp (im(p + ihz - icot), 
where co is the oscillatory frequency, and m and are the azimuthal and longitudinal wavenumbers, 
respectively. The efi'ect of line-tying of the flux tube footpoints in the photosphere is introduced 
by a quantization of the longitudinal wavenumber, 

K = Y' for " = 1,2,3..., (4) 

with L the cylinder length. In this model, n = I corresponds to the fundamental mode. Then, 
Equations Q-© become, 

QVr = 2/VoVy + ^ (k,Br + iB') , (5) 
Bo ^ '' 

v\ I m \ 

Qv^ = -liVoVr + Y [^zB^ - -B, j , (6) 

QBr = -BoKvr, (7) 



^ , , 1 zm 

^25, = -z5o v. + -Vr + — v„ 

1 y f. -r 



(8) 
(9) 



Since we assume the 13 = approximation, we also have = 0. In Equations Q-®, the prime 
denotes a radial derivative, v\ = B^/ y/JIpo is the Alfven speed squared, and Q = oj - mUo/r is the 
so-called Doppler- shifted frequency. Note that Qi = co - mvo within the flux tube and Q.^ = oj in 
the coronal medium, meaning that in both regions the Doppler-shifted frequency is independent of 
r for the given flow profile. Equations dS])-® can be combined into a single equation for the total 
pressure perturbation, = BqB^/iu, namely 



P'4 + -/r^ik'-'^\PT = 0, 



with 



(10) 



(11) 



Note that Equation (ITOI) is equivalent to Equation (2.3) of iBodo et al.l (|1989D if the effect of the 
longitudinal flow is dropped from their expression. The Lagrangian radial displacement, = 
ivJCl, is related to the total pressure perturbation as follows. 



The general solution of Equation (flOl) in the internal medium is 

Pji = AiJ^ (kr) + A2Y,„ {kj) , 



(12) 



(13) 



where Jm and are the Be ssel functions of the first and second kind of order m, respectively 
(lAbramowitz & Stegunlll972h . while Ai and A2 are constants, and k^ is given by Equation (fTTI) . We 
seek for regular solutions at r = 0, so we impose A2 = 0. On the other hand, the general expression 
for the total pressure perturbation in the corona is 



Functions of the first ai 
addition, A3 and A4 are constants, and k^ is defined as 



(14) 



with H^^ and H\'^^ the Hankel functions of the first and second kind of order m, respectively. In 
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The condition for outgoing waves is fulfilled by s etting A4 = and selecting the appropriate branch 
of h such that % (kjoj) > (see details in, e.g.. lwilsonlll98ll : ICallvlll986l) . 



The dispersion relation is obtained by imposing the continuity of both Pj and at the cylinder 
boundary, i.e., r = a. 



(16) 



Equation (fT6l ) is the basic dispersion re lation whose solut ions we will discuss in this investigation, 
and is equivalent to Equation (2.10) of iBodo et al.l (|l989r ) although with a different notation. Note 
that in the absence of flow. Equation (fT6l) reduces to the wel l -known dispersion rel ation of MHD 
waves in a magnetic cylinder (e.g., Edwin & Roberts! 1 19831 : iGoossens et al.l l2009h . In order to 
characterize the solutions of Equation (fT6l) . it is useful to consider their wave properties. In the 
general situation and for fixed and real h, m, and vq, complex values of the frequency, oj = 0)^ + iy, 
are expected. If y < we have a damped mode, whereas y > corresponds to a mode whose 
amplitude gr ows in time , i.e., an unst able mode. Solutions with 7 = w ere called neutrally stabl e 
solutions by iBodo et al. (1989, 1996 ) and correspond to trapped waves fedwin & Roberts 1983 ). 
In such a case, solutions with > are body waves and those with < are surfac e waves 
(which were respectively called reflected and ordinary modes by iBodo et al.lll989lll996() . In the 
general case, y ^ and is complex, meaning that modes have mixed properties since no pure 
body or surface modes are possible. If % {k^^ > 3 expected for a small velocity shear, 

the dominant internal wave character is determined by the sign of % {kf^. Regarding the behavior 
of solutions in the external medium, we must note that k^ is also a complex quantity for complex 
CO, so wave modes have mixed propagating and evanescent properties. For this reason, we use the 
more convenient representation of the external solution in terms of Hankel functio ns instead o f 
the modified Bessel functions usually considered for trapped waves with real o) (see ICallylll986|) . 
It is straight-forward to see from Equation (1151) that when co is real and ojp- < k^v\^, kl is real 
and negative. In such a case, one can directly express the function H^l^ {Ka) in terms of K,n (k^a) 
dAbramowitz & Stegunlll972l) . 



In the following investigation, we restrict ourselves to modes with m 0. This is done 
so because the azimuthal flow has no effect on m = modes in the linear regime. In addition, 
it is unlikely that the development of small spatial-scales at the cylinder boundary obtained by 
Terradas et al.l (|2008bl see their Figure 1) are related to sausage-like (m = 0) solutions but more 
likely to mo des with large m, a nd therefore, with small azimuthal length- scales. Hence, the reader 
is referred to lBodo et al.l(|l989h for a study of the m = case. In addition, we only consider modes 
which are radially fundamental. 
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2.1. Thin Tube Approximation 



To shed more light on the nature of solutions of Equation (1161) . we consider the thin tube (TT) 
approximation, i.e., k^a «; 1, which is realistic in the context of transversely oscillating coronal 
loops. Then, we perform a first order expansion for small arguments of the Bessel and Hankel 
functions in Equation (fT6l ) in the case m 0. By this procedure, we obtain the TT version of the 
dispersion relation, namely 



0. 



(17) 



For a very small velocity shear, one can roughly neglect the terms 4vqQ^ and IvoQiin/ \m\, so the 
effect of the shear flow is only kept in Q.^. In such a case. Equation (fTTl) becomes a second-order 
polynomial in to, whose analytical solution is 

nl/2 



OJ 



Pi +Pe 



-mvo 



■I 



V^,, -l-PeV 



(Pi + Pe) 



PiPe 



.2, .2 



-m V, 



(18) 



(Pi +Pe) 

This solution corresponds to the usual forward (cur > 0) and backward (ojr < 0) body waves 
modified by the shear flow. We see that the KHI takes place for a velocity shear at the cylinder 
boundary that satisfies 



VAi 



> 




(19) 



For a velocity shear larger than the critical value, a damped mode (y < 0) and an unstable mode 
(y > 0) are present. We focus on the unstable solution, whose growth rate in the TT case can be 
easily obtained from Equation (fTSi) . namely 



r 



PiPe 



(Pi + Pe) 



2 2 



(PiVi.+PeVL) 
(Pi + Pe) 



1/2 



(20) 



Note that Equation (|20l) is formally identical to Equation (1) of iTerradas et al.l (|2008b|) . although 
they obtained their expression by assuming incompressible perturbations in a discontinuous inter- 
face. One has to bear in mind that Equations (fT^ and are only approximately valid for a very 
small velocity shear and miss some effects such as, for example, the dependence on the sign of m. 
Therefore, both the actual critical velocity shear for destabilization and the growth rate might be 
significantly different. However, Equations (fT9l ) and (|20l) allow us to predict two important results. 
First, the longitudinal wavenumber has a stabilizing effect since the larger k^, the larger the critical 
velocity shear and the smaller y; and second, modes with large m are more unstable than those with 
small m. This conclusion is equivalent to the well-known general result of the KHI in magnetic 
interfaces (e.g.. IChandrasekhaiill96ll : iDrazin & Reidlll98l|) concerning the stabilizing role of the 
wavenumber component along the magnetic field direction. 
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3. RESULTS 

3.1. Instability Regimes 

First, we study the critical velocity shear that leads to the instability of the solutions of the 
dispersion relation (Equation ([761)). Unless otherwise stated, we consider the next set of param- 
eters: Pi/pe = 3, LI a = 100, and n = I. We restrict ourselves to a velocity shear much smaller 
than the internal Alfven speed, i.e., voa/vAi < 0.1. Equation (fT6l) has been solved for different 
values of the azimuthal wavenumber m. As expected, the degeneracy between forward (o^r > 0) 
and backward (ojr < 0) waves is broken by the presence of flow. When the critical value of the 
velocity shear is reached, the phase velocity of both forward and backward waves coincides and 
the KHI appears. We find that, in the range of considered velocity shears, the first unstable modes 
are those with m = +4, and Figure [T] shows their phase velocity. The behavior of other solutions 
with larger m is similar in all cases, the critical velocity shear being shifted towards smaller values 
as \m\ increases. Note that solutions become leaky, i.e., their phase velocity exceeds the external 
cut-ofi" values, denoted by horizontal dotted lines in Figure [U before the KHI takes place. 

We compare in Figure [2ta) the actual critical velocity shear with that obtained from the ap- 
proximation given by Equation (fT91 ). For small |m|, the actual critical velocity is larger than the 
approximated value. For |m| = 4 there is a factor 2, approximately, between the actual value and 
the approximation. The agreement between both values is substantially better as |m| increases. On 
the other hand. Figure |2b) displays the growth rate of the unstable modes with -20 < m < 20 for 
Voa/vAi = 0.05. The actual growth rate is compared to that predicted by Equation (|20l ) in the TT 
case. We see that the approximated value is slightly larger than the actual one. Again, the differ- 
ence between both of them gets smaller as the azimuthal wavenumber increases. It is worth noting 
that results for positive and negative values of m do not show significant differences regarding both 
the critical velocity and the growth rate. The reason for this behavior is that the considered veloc- 
ity shears are very small, thus the degeneracy between positive and negative values of m is only 
slightly broken. For larger velocity shear, a significant dependence on the sign of m appears. 



3.2. Parametric Study 

In this section we study the dependence of the KHI growth rate, y, on different parameters. 
As stated before, since for the considered shear flow velocities there are no significant differences 
between solutions with positive and negative values of m, hereafter we only show the results cor- 
responding to positive values of m. 

First, we consider pi/pe = 3, Lja = 100, and n = I, and study the dependence of y with 
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Fig. 1 . — Normalized phase velocity versus the velocity shear at the cylinder boundary for Pi/pe = 
3, LI a = 100, and n = I. (a) Solution with m = 4. (b) Solution with m = -A. Dotted lines denote 
the external cut-ofF phase velocities, i.e., iolk^ = ±VAe- 
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Fig. 2. — (a) Critical velocity shear at the tube boundary as a function of m. The symbol A denotes 
the actual value obtained by solving the dispersion relation, whereas the symbol x corresponds 
to the approximation given by Equation (fT9l) . (b) Normalized growth rate as a function of m for 
voa/vAi = 0.05. Again, the symbol A denotes the actual solution and the symbol x corresponds to 
the approximation given by Equation (|20l) . In all cases, Pi/pe = 3, L/a = 100, and n = I have been 
considered. 
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the velocity shear for different values of the azimuthal wavenumber, m. These results are plotted 
in Figure [3ta). As was already commented, we obtain that the larger m, the larger y and the 
smaller the critical velocity shear. Next, we investigate the dependence of y on the longitudinal 
wavenumber, k^. We have fixed m = 6 in the following computations, although equivalent results 
are obtained when other values of m are considered. The value of depends on both n and L 
through Equation dH). Figure Ob) displays the results for different values of n. We see that as the 
value of n increases, i.e., k^ grows, the critical velocity shear is shifted toward larger values and 
the growth rate is reduced. This behavior qualitatively agrees with the analytical predictions of the 
TT case (see Equations (fT9l ) and (|20l) ). and indicates that longitudinal harmonics are more stable 
than the corresponding fundamental mode. On the other hand. Figure (Sjc) shows the results for 
different values of L/a. As L/a increases, k, decreases. We therefore obtain that long loops are 
more unstable than short loops for fixed values of m and n. 

Finally, Figure H] displays the critical velocity shear as a function of the density contrast, 
Pi/Pe- Whereas small values of Pi/pe are realistic in the context of transversely oscillating coronal 
loops, larger values (pi/pe ~ 200) are usually considered to model filament and p rominence fine 



struct i ires (threads). Since k i nk waves have b een detected in filament threads (e.g. JOkamoto et al 



20071 : iTerradas et al.ll2008al : iLin et al.ll2009h . it seems appropriate to extend our present study to 



the context of prominences. We obtain that the critical velocity shear of modes with small m 
increases dramatically as the density contrast grows, although this effect is not so important when 
large values of m are taken into account. This dependence on the density contrast can be seen in 
Equation (fT9l ). These results suggest that the KHI is much more difficult to develop in filament 
thread conditions in comparison to coronal loop conditions, meaning that filament threads might 
be more stable than coronal loops when shear flows are present. 



3.3. Energy considerations 



Results from previous Sections indicate that the larger the azimuthal wavenumber m, the more 
unstable the mode. However, we also have to bear in mind that modes with large m could be very 
difficult to excite by an externa l disturbance, say a fl are, and therefore would weakly contribute 
to the development of the KHI. ITerradas et al.l (|2007h studied how the energy from an initial dis- 
turbance is distributed in different eigen modes of a flux tube. By considering the mathematical 
method by iRuderman & Roberts! (120061) . these authors analytically computed the amplitude and 
the amount of energy de posited in each eigenm ode. Although this calculation depends on the form 
of the initial disturbance. ITerradas et al.l(|2007b concluded that the trapped energy of the eigenmode 
dramatically decreases as the azimuthal wavenumber m is increased. Typically, the difference in 
energy trapped by modes with consecutive m is around three orders of magnitude (see Figure 7(b) 



-12- 




Fig. 3. — Normalized growth rate as a function of the velocity shear at the cylinder boundary, (a) 
Results for m = 4, 5, 6, and 7 with pi/pe = 3, L/o = 100, and n = I. (b) Results for n = 1,2, and 
3 with pi/pe = 3, L/fl = 100, and m = 6. (c) Results for L/a = 50, 100, and 200 with pi/pe = 3, 
m = 6, and n = 1. 
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Fig. 4. — Critical velocity shear at the tube boundary as a function of pi/pe for m = 4, 6, 10, and 
20. In all calculations, L/a = 100 and n = \. 



of iTerradas et al.l l2007|) . Let us assume that the trapped energy by the kink (m = 1) mode is 
El. Then, the energy initially deposited in the eigenmode with azimuthal wavenumber m, namely 
E,„ (t = 0), can be roughly estimated as 



E,n {t = 0)~EiX 10 



-3(|m|-l) 



(21) 



The reader must be aware that Equation (|2TI) is a highly approximated expression since the actual 
deposited energy depends strongly on the form of the initial disturbance among other effects. So, 
the following calculations should be interpreted from a qualitatively point of view. 



On the other hand and according to lBray & LoughheadI (|1974|) . the wave energy density, e, is 
quadratic in the perturbations, thus 



1 

' = 2 



Po(vj+vl) + -(Bj + Bl + Bl) 



(22) 



If a mode is unstable, its related energy density will increase in time by a factor exp (ly^t), where 
7„, is the growth rate for the azimuthal wavenumber m. Therefore, by adding this factor to Equa- 
tion (|2TI) . we get a very approximated estimation for the time-dependent energy of an unstable 
mode with respect to the kink mode energy. 



E,„ (t) = Em (t = 0)x exp (2y,nt) « £i x 10 



-3(|m|-l) 



exp (2y„,0 



(23) 



The quantity E^ (t) /Ei for a fixed time indicates the effective instability of an eigenmode, since 
it takes into account both the initial energy deposited in the mode during its excitation, and its 
subsequent increase due to the instability. A comparison of the value of E^ (0 /Ei corresponding 
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to different modes indicates which of them are the more effectively unstable a the given time- 
scale, t. For the present application, an appropriate time-scale is t = Pk/2, where Pj, is the kink 
mode period. Half a kink mode period is the maximum time for the KHI to develop before shear 
flows change direction due to the flux tube kink motion, and so our assumption that the flow is 
time-independent is not fulfilled. 

Figure [5] displays Em{Pk/2) IE\ for the modes with 4 < m < 20. We see that although the 
growth rate is larger for modes with large m, they cannot acquire enough energy during half a kink 
period to become more relevant than solutions with small m. This qualitative analysis points out 
that only the first two or three unstable modes could effectively contribute to the development of 
a KHI during the kink motion of the flux tube. However, since {Pkl'^) IE\ takes very small 
values, e.g., £■„, {Pk/2) /Ei « 10"^ for the most effectively unstable mode, the KHI cannot remove 
enough energy from the flow to be relevant for the damping of kink oscillations, at least in the 
linear regime. 
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Fig. 5. — Em (Pk/^) /El of the first unstable modes, considering p,/pe = 3, L/a = 100, and n = I. 



3.4. Stabilization by Magnetic Twist 



Results from previous sections indicate that a KHI should develop during kink oscillations of 
coronal loops. However, as these kind of instabilities have not been reported to date, it is likely that 
some mechanism, not considered here, is able to suppress the K HI in coronal loops. From studies 
of the KHI in simple configurations (e.g.. IChandrasekhaiil 1 96 ih . it is known that the inclusion of a 
component of the magnetic field along the flow direction has a stabilizing effect. In our cylindrical 
configuration, this corresponds to an azimuthal component of the magnetic field, so field lines 
would be twisted. The investigation of the effect of magnetic twist on the MHD eigenmodes of 



a flux tube is mathematically complic ated and has b e en broached by a nu mber of authors using 
severa l approximations. For example, iBogdanI (Il984|) . iBennet et alj (Il999|) . and Erdelyi & Fedun 
(120061) st udied wave propagation in tw isted tubes by tak i ng the incompressible approximation into 



account. Erdelyi & Carted (|2006|) and ICarter & Erdelyil (|2007l |2008|) investigated w ave modes in 
cylind ers with the magnetic twist being constrained within an annulus. whereas £rdelvi & Fedun 
(|2007h only considered m = modes in fully twisted tubes. Finally, iRudermanI (|2007h studied 
non-axisymmetric oscillations of twisted tubes in the thin tube approximation. It is beyond the 
purpose and scope of the present work to include magnetic twist in our cylindrical configuration, 
since it would complicate matters, e.g., by introducing resonances in the system. Instead, we 
consider in this section the more simple Cartesian (slab) geometry, and study how a component of 
the magnetic field along the flow direction affects the KHI. An analogy between the slab case and 
the cylinder case can be performed. 

Let us consider a magnetic slab with half-width a in the .^-direction embedded in an unlimited 
environment. Both the y- and z-directions are unbounded. All symbols in the following equations 
have the same meaning as in previous sections. We assume a mass flow in the y-direction, Uo = 
Uoey. The following density, po, and flow velocity, Uq, profiles are assumed, 

if \x\ < a, ( VqU, if \x\ < a, 



Po = 



Pe, 



if \x\ > a. 



Ua = 



0, if \x\ > a. 



(24) 



The magnetic field is 



Bn = 



(25) 



5o cos or gj, + 5o sin or e^, if |x| < a, 
B^e^, if \x\ > a. 

Thus, for a = njl the internal magnetic field is perpendicular to the flow, whereas for a = the 
internal magnetic field and the flow are in the same direction. We now apply Equations ©-([S]) 
to this configuration and assume perturbations proportional to exp {ikyj + ik^z - icotj, with k^ given 
by Equation & and ky = ml a. This forrn of ky has been assumed in a number of works (e.g., 
HoUweg & YangI 1 19881 : lArregui et al.l 120071) in order to ext end the results of the slab models to 
the cylindrical geometry. In particular, lArregui et al.l (|2007|) showed that the solutions of the slab 
model tend to those of the cylindrical configuration when larger values of ky in comparison with k^ 
are taken into account. In the present case, ky » h as long as Lja » nn/m. 

The general d i spersi on relation of MHD modes supported by the slab model is given in Equa- 
tion (7) of ljoardeil toO^ . In our present notation and for jS = 0, the dispersion relation becomes. 



k^vli sin^ (e + a)] 



tanh {%a) 
coth {%a) 



+ 



p,%(co^-klvl^) = 0. 



(26) 



with 
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where Cli = co - mvq, = ky + k^, and 6 = arctan (ky/k^. The tanh term in Equation (|26|) 
corresponds to symmetric solutions whereas the coth term stands for anti-symmetric solutions 
with respect to the x-axis. In the case a = n/2 and vq = 0, Equation (|26l) reduces to the dispersion 
relation derived by lRae & Roberts! (|1983|) . 



Figure |6ta) shows a comparison of the previously obtained growth rates in the cylinder case 
with those computed in the slab case by solving Equation (|26l) . We see that the growth rates 
are slightly larger in the slab case than in the cylinder case but both solutions are in reasonable 
ag reement, the differen ce being reduced as m increases. This behavior coincides with the results 
of lArregui et al.l (|2007[) for the real part of the frequency. On the other hand, symmetric and anti- 
symmetric modes in the slab case have very similar growth rates, so hereafter we restrict ourselves 
to symmetric modes for simplicity. Next, Figure [6tb) displays the growth rate of symmetric modes 
as a function of a for voa/vAi = 0.1 . We see that, depending on m, there is a critical a which 
suppresses the KHI. The critical a clusters towards a limit value as m increases. To obtain an 
analytical expression for this limit value of a, let us consider the situation m » ha. In such a case, 
the length-scale in the y-direction is much smaller than those in the x- and z-directions, and the 
solutions of the slab configuration can be approximated by those of a single interface. In addition, 
if we assume incompressible pertu rbations, the dispersion relation is (|joarder & Satya Narayanan 
2000l : ljoarder & Nakariakovlbooel) . 



Pi _ k\l. sin^ (6 + a)) + p, (co^ - k^^^) = 0, 
whose exact analytical solution is 



to = 



Pi +Pe 



-mvo + 



(piv2.F sm\e + a) + p.vlk^) 



Pipe 



2, ,2 



(Pi + Pe) 



rtn V, 



1/2 



(28) 



(29) 



(Pi +Pe) 

Note the remarkable resemblance between Equations (fTSl) and (|29l ). although the former was ob- 
tained for a cylinder in the TT and /3 = approximations, and the latter corresponds to a magnetic 
interface in the incompressible limit. In fact, it is possible to directly obtain Equation (l29l) from 
Equation (fTSl) by simply adding the contribution of an azimuthal component of the magnetic field 
considering its value at the cylinder edge. Nevertheless, we have derived Equation (|29l ) from the 
expressions of the Cartesian case for the sake of a more rigorous method. 



Finally, it is straight-forward to derive from Equation (l29l) the condition for stability, 

v2 ,9 9 



sin (O + a) = 



m 



(Pi+Pe) (m^ + kja^) 



voa 

VAi 



ky 



(30) 



By performing the limit m » fc^a one obtains that the critical a satisfies. 



cos a 



Voa 



iPi + Pe) \ VAi 



(31) 
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Fig. 6. — (a) Normalized growth rate as a function of m for voa/vAi = 0.05. The symbol A denotes 
the solution of the cylinder case, whereas symbols x and O correspond to the symmetric and anti- 
symmetric slab solutions, respectively, (b) Normalized growth rate of symmetric slab solutions as 
a function of a for m = 4, 6, 10, and 20 with voa/vAi = 0.1. The vertical dotted line denotes the 
critical angle for m ^ oo given by Equation (|3TI) . In all cases, pjpe = 3, L/a = 100, and n = I 
have been considered. 
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The value obtained from Equation (|3TI) is indicated in Figure [6tb) by a vertical line, showing a 
good agreement with the results. For a smaller than the limit critical value, the KHI is completely 
suppressed for any value of m. For voa/vAi = 0.1, the limit critical angle is a a: 87deg, which 
corresponds to a small inclination of the magnetic field with respect to the z-axis, i.e., a very small 
azimuthal component of the magnetic field in the analog cylindrical case. So, a very "weak" twist 
of magnetic field lines can be enough to prevent the triggering of the KHI in a cylindrical coronal 
loop, the azimuthal component of the magnetic field needed to suppress the KHI being probably 
much smaller than the critical value of the pinch instability. A weak twist of magnetic field lines 
is very likely and realistic in the context of coronal loops. Although further studies are needed, 
magnetic twist seems to be a consistent explanation for the absence of the KHI in the observations 
of transversely oscillating loops. 



4. CONCLUSION 

In this paper, we have studied the linear KHI in straight, untwisted, and cylindrical coronal 
magnetic flux tubes due to the presence of azimuthal shear flows at the tube boundary. We have 
found that for realistic values of the velocity shear, fluting-like MHD modes with large azimuthal 
wavenumber and small longitudinal wavenumber are unstable. A qualitative analysis based on 
energy considerations suggests that only the first unstable fluting modes can develop a KHI during 
time-scales related to kink motions of the flux tube. Since these kind of instabilities have not been 
reported by observers to date, we considered an analog in Cartesian geometry and obtained that 
the twist of magnetic field lines could completely suppress the KHI in twisted coronal flux tubes. 
The magnetic twist that can suppress the KHI is very weak and probably much smaller than the 
critical twist of the pinch instability, meaning that the magnetic stability of the flux tube would not 
be compromised. Further investigations on this issue are needed. 

Some efi'ects not considered here could be studied in future works. The next logical step is 
to assume a time and azimuthal dependence of the flow in order to describe more realistically the 
actual motions of a kink oscillation. This investigation should be numerically performed. Other 
effects are, for example, to consider the jS ^ case, or to take into account the presence of a 
transitional layer between the loop and the corona. 
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